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Abstract 



We present the results of a lattice QCD calculation of the average up-down and 
strange quark masses and of the light meson pseudoscalar decay constants with 
Nf = 2 dynamical fermions. The simulation is carried out at a single value of the 
lattice spacing with the twisted mass fermionic action at maximal twist, which guar- 
antees automatic C(a)-improvement of the physical quantities. Quark masses are 
renormalized by implementing the non perturbative RI-MOM renormalization pro- 
cedure. Our results for the light quark masses are m^f(2 GeV) = 3.85 ± 0.12 ± 0.40 
MeV, mf (2 GeV) = 105±3±9 MeV and m s /m ud = 27.3±0.3±1.2. We also obtain 
fx = 161.7±1.2±3.1 MeV and the ratio fx / f-w = 1.227±0.009±0.024. From this ra- 
tio, by using the experimental determination of T(K — > fiv^)) /T(tt — > ^1/^(7)) and 
the average value of \V u d\ from nuclear beta decays, we obtain \V US \ = 0.2192(5)(45), 
in agreement with the determination from K\$ decays and the unitarity constraint. 

1 Introduction 

In this paper we extend to the kaon sector our previous lattice study of the pion mass 
and decay constant pQ. We present a determination of the light quark masses, strange 
quark mass m s and the average up-down quark mass m u d, of the kaon pseudoscalar decay 
constant fx, and of the ratio fx/ fn- We have simulated the theory with Nf = 2 dynamical 
quarks, taken to be degenerate in mass, and two valence quarks. In order to investigate 
the properties of the K meson, we consider in the present analysis a partially quenched 
setup, namely we take the valence quark masses fii and H2 different in value between each 
other and different from the sea quark mass fi$- 

The strategy of the calculation is the following. We first compute the pseudoscalar 
meson masses and decay constants for different values of the sea and valence quark masses, 
and study their mass dependence. We then use the experimental values of the ratios M n / f n 
and Mk/Mk to determine the average up-down and the strange quark mass respectively. 
The lattice spacing is fixed from The results obtained for the quark masses are finally 
used to evaluate fx and the ratio fx/ fv 

The calculation is based on a set of gauge field configurations generated with the tree- 
level improved Symanzik gauge action at f3 — 3.9, corresponding to a = 0.087(1) fm 
(a -1 ~ 2.3 GeV) |1], and the twisted mass fermionic action at maximal twist. We have 
simulated 5 values of the bare sea quark mass, 

an s = {0.0040, 0.0064, 0.0085, 0.0100, 0.0150} , (1) 

and computed quark propagators for 8 values of the valence quark mass, 

a//i, 2 = {0.0040, 0.0064, 0.0085, 0.0100, 0.0150, 0.0220, 0.0270, 0.0320} . (2) 

The first five masses are equal to the sea quark masses, and lie in the range 1/6 m s ^ 
A*i,2 ~ 2/3 m s , where m s is the physical strange quark mass, while the heaviest three are 
around the strange quark mass. 
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We implement non-degenerate valence quarks in the twisted mass formulation of lat- 
tice QCD as discussed for instance in refs. [2J [3]. We introduce two twisted doublets of 
degenerate valence quarks, (u±, di) and («2, c^), with masses fii and \Li respectively, and 
simulate charged mesons u\di and d\Ui- Within each doublet, the two valence quarks are 
regularized in the physical basis with Wilson parameters of opposite values {r u = — = 1). 

At each value of the sea quark mass we have computed the two-point correlation func- 
tions of charged pseudoscalar mesons, with both degenerate and non degenerate valence 
quarks, on a set of 240 independent gauge field configurations, separated by 20 HMC trajec- 
tories one from the other (each trajectory being of length 1/2). To improve the statistical 
accuracy, we have evaluated the meson correlators using a stochastic method to include all 
spatial sources. The method involves a real stochastic source (Z(2)-noise) for all colour and 
spatial indices at one Euclidean time slice randomly moved when passing from one gauge 
configuration to another. This "one-end" method is similar to that pioneered in ref. [4J 
and implemented in ref. [5]. Statistical errors on the meson masses and decay constants 
are evaluated using the jackknife procedure, by decimating 10 configurations out of 240 in 
each jackknife bin. Statistical errors on the fit results, which are based on data obtained 
at different sea quark masses, are evaluated using a bootstrap procedure. Further details 
on the numerical simulation can be found in [6]. 

The use of twisted mass fermions in the present calculation turns out to be beneficial 
in several aspects 012]: i) the pseudoscalar meson masses and decay constants, which 
represent the basic ingredients of the calculation, are automatically improved at 0(a); ii) 
once maximal twist is realized, the physical quark mass is directly related to the twisted 
mass parameter of the action, and it is subject only to multiplicative renormalization; iii) 
the determination of the pseudoscalar decay constant does not require the introduction of 
any renormalization constant, and it is based on the relation 

fps = (mi + M2) — — • ( 3 ) 

Concerning the size of discretization effects, it is worth noting that, since the two valence 
quarks are regularized in the physical basis with Wilson parameters of opposite values, 
the meson mass Mp S differs from its continuum counterpart only by terms of 0(a 2 ^) and 
0(a 4 ), whereas fps differs from its continuum limit by terms of 0(a 2 ) [8j [9]. Therefore, 
at 0(a 2 ) the cutoff effects on Mp S and fpg are as in a chiral invariant lattice formulation. 

The meson mass Mp$ and the matrix element |(0|P 1 (0)|P)| have been extracted from 
a fit of the two-point pseudoscalar correlation function in the time interval t/a £ [10,21]. 
In order to illustrate the quality of the data, we show in fig. [I] the effective masses of 
pseudoscalar mesons, as a function of the time, in the degenerate cases — A*i = A*2- 

1 Strictly speaking, automatic O(a) improvement was proved in [71 [5] to hold in a unitary as well as in 
a mixed action framework. Actually the same proof goes through also in the present partially quenched 
setup. The reason is that all the symmetries entering the discussion of the renormalizability and O(a) 
improvement are valid for generic values of the masses of the various valence and sea quarks. 
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2 Quark mass dependence of pseudoscalar meson mas- 
ses and decay constants 

The determination of the physical properties of K mesons requires a study of the quark 
mass dependence of the corresponding observables over a large range of masses, extending 
from the physical strange quark down to the light up and down quarks. In this work, we 
study the quark mass dependence of the pseudoscalar meson masses and decay constants 
by investigating two different functional forms. The first one is the dependence predicted 
by continuum partially quenched chiral perturbation theory (PQChPT), whereas in the 
second case we consider a simple polynomial dependence. For a recent precision study of 
the quark mass dependence of meson masses and decay constants in the partially quenched 
theory with Nf = 2 dynamical fermions see also ref. [10J. 



2.1 PQChPT fits 



Within PQChPT we consider the full next-to-leading order (NLO) expressions with the 
addition of the local NNLO contributions, i.e. terms quadratic in the quark masses, which 
turn out to be needed for a good description of the meson masses and decay constants up 
to the region of quark masses around the strange quark mass. The PQChPT predictions 
have been derived in ref. [TT] and can be written in the form 



1 + 6 (6~6) In 2fr _ 6(6~6) In 26 | 



(6-6 



(6-6 



+a v 62 + as £s + a vv + a ss £§ + a V s 62 £s + a VD £,012 



(4) 



3 



fps((*S, 1*1,1*2) = f 



1 - 6s In 2d, - 6s In 26s + %fff In (|) + 



+{b v + 1/2) 62 + (6s - 1/2) £s + b vv & + &ss Cs + &vs 62 £s + b VD £ 2 12 

where 6 = 2B fj li / '{Air f f , 6j = B (^ + /i i )/(4 7 r/) 2 and £ Dij = B (^ - /i i )/(4?r/) 2 . The 
parameters £?o and / are the low energy constants (LECs) entering the chiral Lagrangian 
at the LO0, whereas ay, as, by and bs are related to the NLO LECs pj] by 



ay = 4«g — 2«£ 



05 = 8«6 — 4a 4 



«5 



2«m 



(5) 



The quadratic terms in the quark masses in eq. (jlj) represent the local NNLO contributions. 
The corresponding chiral logarithms at two loops in the partially quenched theory are also 
known [T2]. They involve, however, a larger number of NLO LECs whose values, in the 
Nf = 2 theory, cannot be fixed from phenomenology. Introducing their contribution in the 
fit would increase significantly the number of free parameters, thus limiting, at the same 
time, the predictive power of the calculation. 

In the limit of degenerate valence quark masses, /xi = /x 2 = f*v, eq. (jlj) is finite and 
reduces to 



Mp S {fis, Hv, liv) = 2 B fi v • [1 + (2 £y - &) In 2£y + (ay + 1) £y + (a s - 1) 6 

+ a VV £y + a ss £s + a vs £v 6 
fps(p-S, t*v, l*v) = f ■ 1-2 £ys In 2f vs + fry £y + &s £s + £y + b ss £s + Vs £v 6 



(6) 



2.1.1 Finite volume corrections 

In a lattice QCD calculation aiming at a percent precision on the physical predictions, the 
impact of finite size corrections cannot be neglected. The lattice in our simulation has 
spatial extension L = 24a ~ 2.1 fm, and the pseudoscalar meson mass at the lightest value 
of the quark mass is such that MpsL ~ 3.2. Since we have not performed yet a systematic 
study of non-degenerate meson masses and decay constants on different lattice volumes, 
we will estimate the finite size effects by including in the fits the corrections predicted by 
one-loop chiral perturbation theory, which, in the partially quenched expressed 
by pj 1 



Mp S (/j,s, m,IMi; L) = Mp S (fi S , Hi,fJ>2) 
fps(fis,Vi>i,H2]L) = fpsifis, 1*1,1*2) ■ 



1 + 



6 (£s - 6) 9i(L, 6) £2 (£s - 6) 9i{L, 6) 



(6-6 



(6-6 



1 - 6s 9i (L, 6s) - 6s 9i (L, 6s) + 



62 - £s 

2(6-6: 



(6 91 (£,6) -6£i(£,6)) + 



+^ (6? - 6) 6) + J (6? - 6) 6) 



(7) 



2 The pseudoscalar decay constant / is normalised such that f„ = 130.7 MeV at the physical pion mass. 
3 We thank D.Becirevic and G.Villadoro for having provided us with the expression of finite volume 



corrections to Mp S (fXs, Mi> M2) which is not given in ref. 
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The functions g s (s — 1,2) in eq. ([7]) are defined as 

(4tt) 3 / 2 



9s(L,M 2 



r(s- 1/2) £ S _ 1/2 (L,M 2 ) 



(M 2 ) 2 - s 

where M is the pseudoscalar meson mass at the LO, M 2 = 2_B At = (47r/) 2 £, 
and $(r) is the elliptic theta function 

oo 



(9) 



(10) 



The limits of eq. ((Tj) in the case of degenerate valence quark masses, fix = // 2 = A*y, can 
be obtained by using the identity 



M 



c/M 2 



g s (L, M 2 ) = -(2 - s) ~g s (L, M 2 ) - ~g s+1 (L, M 2 ) 



and are given by 



fps(jis,fiv,fiv',L) = f PS (ns,Vv,Hv) ■ [1 - 2£ vs g 1 (L,£vs)] ■ (12) 
2.2 Polynomial fits 

The inclusion of the local NNLO contributions in the PQChPT predictions expressed by 
eq. (j3J) is required by the observation that the pure NLO predictions are not accurate 
enough to describe the quark mass dependence of pseudoscalar meson masses and decay 
constants up to the region of the strange quark. However, not having considered the full 
NNLO chiral predictions, we regard eq. (|3J) mostly as an effective description of the quark 
mass dependence of these observables. In order to evaluate the associated systematic 
uncertainty, we also consider in the analysis an alternative description based on a simple 
polynomial dependence on the quark masses, for both the pseudoscalar meson masses and 
decay constants: 



' 1 + a V £l2 + a S £,S + a VV £12 + a SS d + a VS £l2 £s + a VD ^D12 j 

fps{vs,lH,tki) = /• [l + (bv + 1/2) £12 + (6s- 1/2) & + b vv e i2 + bss 

+bvs 62 £s + ^Di2 



(13) 



Note that, though we are adopting in eq. (fl3|) the same notation for the coefficients of the 
chiral expansions as in eq. (j3J), the physical meaning of these coefficients, i.e. their relation 
to the derivatives of Mp S and fps with respect to the quark masses, is actually different. 
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It also worth observing that, in the case of the polynomial fits f fl3|) . a change in the values 
of the LECs / and B only amounts to a redefinition of the fit parameters of Mp S and fpg 
respectively. Therefore, in this case, the two fits are independent one from the other. The 
differences between the results obtained by performing either chiral or polynomial fits will 
be included in the final estimates of the systematic errors. 

3 Chiral extrapolations 

The input data in the present analysis are the lattice results for the pseudoscalar meson 
masses and decay constants obtained at each value of the sea quark mass, with both 
degenerate and non degenerate valence quarks. We exclude from the fits the heaviest 
mesons having both the valence quark masses in the strange mass region, namely with 
^1,2 = {0.0220, 0.0270, 0.0320}. Overall, we have considered therefore 150 combinations 
of quark masses for both the meson masses and the decay constants. The full sets of 
results are collected in tables H] and [5] of the appendix. The number of free parameters in 
the combined fit of Mp S and fpg is 14, but a first analysis shows that some of them, in the 
various cases, are compatible with zero within one standard deviation, and are kept fixed 
to zero in the final estimates of the fit parameters (see table [1]) . 

In order to extrapolate the pseudoscalar meson masses and decay constants to the 
points corresponding to the physical pion and kaon, we have considered three different fits: 

• Polynomial fit: a polynomial dependence on the quark masses is assumed for the 
pseudoscalar meson masses and decay constants, according to eq. ([TBI . 

• PQChPT fit: the pseudoscalar meson masses and decay constants are fitted ac- 
cording to the predictions of PQChPT expressed by eq. (jl]) to which we add the 
finite volume corrections of eq. (J7j). 

• Constrained PQChPT fit : this fit, denoted as C-PQChPT in the following, de- 
serves a more detailed explanation. The main uncertainty in using eqs. (jlj) and ( {TBI 
to effectively describe the quark mass dependence of Mp S and fps is related to the 
extrapolation toward the physical up and down quark masses. On the other hand, 
we have shown in ref. [1] that pure NLO ChPT, with the inclusion of finite volume 
corrections, is sufficiently accurate in describing the lattice results for both the pseu- 
doscalar meson masses and decay constants when the analysis is restricted to our 
lightest four quark masses in the unitary setup (i.e. Hi = /i 2 = /ig). In order to take 
advantage of this information, when performing the C-PQChPT fit we first determine 
the LO parameters B and / and the NLO combinations ay + a$ and by + b$ from a 
fit based on pure NLO ChPT performed on the lightest four unitary points. In other 
words, we repeat here as a preliminary step the same analysis done in ref. [JJ, but 
on the smaller statistical sample of data used for the present studyS In this way we 

4 Note that in the limit = ji2 = Ms, an d when all the coefficients of the quadratic terms are sent to 
zero, the PQChPT expressions (f4]), as well as the hnite volume corrections expressed by eq. reduce to 
the pure NLO ChPT predictions used in the chiral fit of ref. pQ. 
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determine 



2aB = 4.82(10) , af = 0.0552(12) , 

a v + a s = 0.80(23) , b v + b s = 0.62(24) . (14) 

These results, are perfectly consistent, at the level of ~ 1.5 a, with those obtained 
in ref. pp. By using the constraints of eq. (JT41 . the other parameters entering the 
chiral expansions of Mp S and fps are then obtained from a fit to eq. (jl]) over the non 
unitary points. For consistency with the previous unitary fit, we exclude also in this 
case from the analysis the data at the highest value of sea quark mass, afis = 0.0150. 

In table [1] ( "All data" ) we collect the results obtained for the fit parameters in the three 
cases: polynomial, PQChPT and C-PQChPT fits. In the last line we also quote the 
corresponding values of the x 2 P er degree of freedom. From these values we see that, though 
the quality of the fit is better in the polynomial case, all three analyses provide a good 
description of the lattice data, in the whole region of masses explored in the simulation. 
This is only true, however, if the terms quadratic in the quark masses are taken into 
account. 

A potential problem in the partially quenched theory is the divergence of the chiral log- 
arithms in the limit in which the light valence quark mass goes to zero at fixed sea quark 
mass (see eq. fll])). This divergence does not affect the extrapolation of the lattice results to 
the physical point, since the sea and the light valence quark masses are degenerate in this 
case. However, in order to verify that this unphysical behaviour of the partially quenched 
chiral logarithms does not modify the result of the extrapolation, we have repeated the 
analysis by restricting both the polynomial and the chiral fits to the 30 quark mass combi- 
nations (26 in the case of the C-PQChPT fit) that, satisfying the constraint \ii > /ii = fis, 
are not affected by the dangerous chiral logarithms. The results obtained for the free pa- 
rameters of these fits are also shown in table [1] (last three columns) . By comparing these 
results with those obtained by using the full set of data, we find some differences in the 
estimates of the coefficients of the quadratic terms, particularly those involving the sea 
quark mass (ays, ass, • • •)• These differences reflect the relative importance in the fit of 
the various quadratic terms in the different quark mass regions. For instance, in the case 
of the highest sea quark mass, afis = 0.0150, only 4 out of 30 combinations of masses are 
included in the fit restricted by the condition ji 2 > Hi = Hs- On the other hand, when we 
compare the results for the extrapolated physical quantities (am ud , am s , af n , . . .) obtained 
from the two fits, we find that they are almost indistinguishable (see table [2]). This is 
reassuring, as it shows that the effects of potentially divergent chiral logarithms are well 
under control in our analysis. 

The mass dependence of the pseudoscalar meson masses and decay constants is illus- 
trated in fig. [21 where we also compare the lattice data with the results of the polynomial, 
PQChPT and C-PQChPT fits. We have shown in the plots the cases in which one of the 
valence quark mass (/ii) is equal to the sea quark mass, and the results are presented as a 
function of the second valence quark mass {^2)- The points corresponding to the physical 
pion and kaon are thus obtained by extrapolating/interpolating the results shown in fig. [2] 
to the limits [i\ — * m u d and [ii — > m s . 
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Figure 2: Lattice results for a 2 Mp S (top), a 2 Mp S /^(afii + a/z 2 ) (center) andafps (bottom) 
as a function of the valence quark mass afj, 2 , with a/ii = a^s- The solid, dashed and dotted 
curves represent the results of the polynomial, PQChPT and C-PQChPT fits respectively. 
For better clarity, results at only two values of the see quark mass have been shown in the 
center plot. 
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All data 


Only ji 2 > [i\ = 


Us 


r re 


Polynomial 






Polynomial 




p pnpv.PT 


2aB 


4.59(3) 


4.79(6) 


4.82(10) 


4.55(6) 


4.86(12) 


4.82(10) 


af 


0.0607(6) 


0.0577(6) 


0.0552(12) 


0.0606(9) 


0.0574(14) 


0.0552(12) 


a v 


-0.63(7) 


2.37(10) 


2.15(18) 


-0.52(16) 


1.91(15) 


2.15(18) 


a s 


0.0 


-1.44(10) 


-1.35(12) 


0.0 


-1.04(37) 


-1.35(12) 


7 

K 


2.66(4) 


0.68(5) 


0.86(8) 


2.56(13) 


0.49(12) 


0.75(8) 


bs 


0.86(13) 


-1.22(15) 


-0.25(23) 


1.03(15) 


-0.94(34) 


f\ "1 O / C\ A \ 

-0.13(24) 


avv 


2.6(2) 


-9.3(3) 


-8.3(6) 


2.3(5) 


-7.8(18) 


-5.8(7) 


ays 


0.0 


7.6(4) 


6.9(3) 


0.0 


6.0(38) 


0.0 


ass 


0.0 


0.0 


0.0 


0.0 


0.0 


5.9(7) 


CtVD 


-0.6(1) 


-3.8(2) 


-3.2(3) 


-0.9(6) 


-2.6(21) 


-5.1(4) 


7 

bvv 


-4.0(2) 


1.2(2) 


0.9(1) 


-4.1(8) 


0.0 


2.3(5) 


bvs 


0.0 


6.0(6) 


3.7(12) 


0.0 


7.1(21) 


0.0 


bss 


0.0 


0.0 


-5.3(14) 


0.0 


0.0 


-2.0(6) 


bvD 


-3.7(2) 


-3.8(2) 


-3.0(3) 


-2.6(6) 


0.0 


-3.1(6) 


X 2 /d.o.f. 


0.38 


1.34 


1.11 


0.28 


0.40 


0.78 



Table 1: Values of the fit parameters as obtained from the polynomial, PQChPT and 
C-PQChPT fits (see text for details), by analysing all combinations of quark masses or only 
the combinations satisfying the constraint Hi>\L\ = fis- In the last line, the corresponding 
X 2 per degree of freedom are also given. 



In order to illustrate the impact of finite volume corrections in the PQChPT fits, we 
compare in fig. [3] the best fit curves for the pseudoscalar meson masses and decay constants 
as obtained with or without including these corrections. In the plots the differences between 
the two curves are barely visible. Obviously, a different question is whether the theoretical 
formulae based on ChPT can accurately describe at the NLO the dependence of Mp S and 
fps on the lattice volume. We postpone this issue to a future investigation, in which 
we plan to better quantify the systematic error due to finite size effects by extending 
the calculation of light pseudoscalar meson masses and decay constants on lattices with 
different spatial sizes. 

By having determined the fit parameters, we are now ready to extrapolate eqs. (j3j) and 
( TT3"j) to the physical pion and kaon. We follow the procedure outlined in sect. [TJ we use the 
experimental values of the ratios M n /f n and Mx/M n to determine the average up-down 
and the strange quark mass respectively. Once these masses have been determined, we use 
again eqs. (jlj) and f[13l) to compute the values of the pion and kaon decay constants as well 
as their ratio fx/ f^u 

5 In order to account for the electromagnetic isospin breaking effects which are not introduced in the 



9 



PQChPT Fit 




PQChPT Fit 



0,09 - 



0,08 



0.07 - 



0,06 



0.1)00 




an s =an,=0.0040 
an s =an 1 =0.0064 
aji s =a|Xj=0.0085 
an s =an.j=0.0100 
an s =aM,,=0.0150 

— With FSE corr. 

— Without FSE corr. 



0,005 0,010 0,015 



0,020 

am 



0,025 



0,030 



0,035 



0,040 



Figure 3: PQChPT fits of the pseudoscalar meson masses and decay constants performed with 
(solid lines) and without (dashed lines) including the finite volume corrections of eq. (0). The 
results are shown as a function of the valence quark mass a^2, with apL\ = afis- 



In table [2] we collect the values of the quark masses, meson masses and decay constants, 
in lattice units, as obtained from the three fits by analysing all combinations of quark 
masses or only the combinations that satisfy the constraint /i2 > A*i = ^s- Note that the 

lattice simulation, we use as "experimental" values of the pion and kaon mass the combinations |14j 
(Ml) QCD = Mlo , (M 2 K ) QCD = i [M|c + M| + - (1 + A E )(M* + - M*,)] 

with Ae = 1. 
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All data 


Only fi 2 > Hi = 




Fit 


Polynomial 


PQChPT 


C-PQChPT 


Polynomial 


PQChPT 


C-PQChPT 


am ud ■ 10 3 


0.90(2) 


0.86(2) 


0.79(4) 


0.91(3) 


0.84(5) 


0.79(4) 


am s 


0.0243(5) 


0.0235(5) 


0.0218(10) 


0.0243(7) 


0.0234(12) 


0.0217(10) 


m s /m ud 


26.9(1) 


27.4(2) 


27.5(3) 


26.7(2) 


27.9(2) 


27.4(3) 


aM w 


0.0642(6) 


0.0632(6) 


0.0610(12) 


0.0642(9) 


0.0629(14) 


0.0610(12) 


aM K 


0.235(2) 


0.232(2) 


0.224(4) 


0.235(3) 


0.231(5) 


0.224(4) 


aU 


0.0622(6) 


0.0612(6) 


0.0591(11) 


0.0622(8) 


0.0609(13) 


0.0591(11) 




0.0756(7) 


0.0744(7) 


0.0730(11) 


0.0755(8) 


0.0747(11) 


0.0731(12) 


fjtj fn 


1.216(3) 


1.215(4) 


1.236(8) 


1.214(8) 


1.225(11) 


1.238(7) 



Table 2: Values of the quark masses, meson masses and decay constants in lattice units as 
obtained from the polynomial, PQChPT and C-PQChPT fits by analysing all combinations 
of quark masses or only the combinations satisfying the constraint /i 2 > fii = Hs- 



values of the quark mass ratio m s /m u d and of the ratio of decay constants fx/fn, being 
dimensionless and well normalised quantities, are obtained at this step without need of 
fixing the scale nor of introducing the quark mass renormalization constant. For these 
quantities, therefore, the results presented in table [2] already represent physical predictions 
of the calculation. 

As a further investigation, we have studied how the results for the quark masses and 
decay constants change when the analysis is performed only on mesons with degenerate 
valence quarks. In this case, we find values of quark mass in good agreement with those 
given in table (T5]), whereas for fx and fx/ fn we obtain results that are larger by about 
5% than those quoted in the table. This reflects the fact that the mass difference between 
valence quarks represents only a small effect in meson masses, while it turns out to be 
relevant in decay constants, at the present level of accuracy, as shown by the contribution 
of the a V D and byD terms respectively in the simple polynomial fits (see tabled]). 

4 Physical results 

In order to convert into physical units the results obtained for the strange quark mass and 
the kaon decay constants we fix the scale within each analysis (polynomial, PQChPT and 
C-PQChPT fits) by using f n as physical input. This choice deserves some discussion. By 
looking at table we see that the value of the pion decay constant in lattice units as 
obtained from the C-PQChPT fit is in agreement, at the level of 1.4 a, with the result 
of our previous study, af n = 0.0576(7) pQ. Indeed, from the present analysis we obtain 
the estimate a = 0.089(2) fm, to be compared with the determination a = 0.087(1) fm of 
ref. pQ. We also find that the estimate of the lattice spacing obtained from the C-PQChPT 
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analysis coincides with the one derived from the pure NLO ChPT analysis performed over 
the lightest four unitary points. This is expected, since as explained before the NLO 
unitary fit over the four lightest quark masses is used as a constraint in the C-PQChPT 
analysis, and the effect of the quadratic terms which are left out in the first fit is negligible 
in the evaluation of /„-. We then conclude that the difference between the determination 
a = 0.089(2) fm and the one given in ref. [1] is a purely statistical effect and, as such, 
is properly accounted for by the quoted statistical errors. In the analyses based on the 
PQChPT and polynomial fits, instead, we obtain the estimates a = 0.092(2) fm and 
a = 0.094(1) fm respectively. In this case, the differences with respect to the C-PQChPT 
determination, which are at the level of 3% and 6% respectively, have a systematic origin 
related to the uncertainty in the chiral extrapolation. 

As mentioned before, rather than choosing a common estimate of the scale for the 
polynomial, PQChPT and C-PQChPT analyses, we prefer to fix the scale by relying on 
the determination of af n as obtained within each separate fit. This choice has the impor- 
tant advantage that also the pion and kaon masses are fixed in this way to their physical 
values within each fit, since the experimental results for the ratios M n /f n and Mx/M n 
have been used to determine the light and strange quark masses. Therefore, the absolute 
normalization of the fit functions describing the quark mass dependence of both the meson 
masses and the decay constants is always correct, independently of the assumptions done 
on the chiral behaviour. As a result, we find that the systematic differences among the 
various determinations of am s and a fx given in table [2} which are at the level of 6% and 
2% respectively, reduce by approximately a factor of two when the results are converted 
in physical units. Nevertheless, in the case of the polynomial and PQChPT fits we conser- 
vatively add in the calculation of the dimensionful quantities a 6% and 3% of systematic 
error coming from the different estimates of the scale. 

The determination of the physical strange and up-down quark masses also requires 
implementing a renormalization procedure. The relation between the bare twisted mass at 
maximal twist, fi q , and the renormalized quark mass, m q , is given by 

fn q {fi R ) = Z m (g 2 , a/i R ) fi q {a) , (15) 

where fip is the renormalization scale, conventionally fixed to 2 GeV for the light quarks. 
Z m is the inverse of the flavour non-singlet pseudoscalar density renormalization constant, 
Z m = Zp . We have used the (9(a)-improved non-perturbative RI-MOM determination of 
Z P , which gives Zf J - MOM (l/a) = 0.39(1)(2) at (3 = 3.9 [15], and converted the result to the 
MS scheme at the scale [ir = 2 GeV by using renormalization group improved continuum 
perturbation theory at the N 3 LO [16J. 

In table [3] we collect the results for the light quark masses and pseudoscalar decay 
constants, in physical units, as obtained from the polynomial, PQChPT and C-PQChPT 
fits. For completeness, we also show in the table the results for the ratios m s /m u d and 
fx I f-K already given in table EJ To be conservative, we only consider from now on the 
results obtained from the analysis of the quark mass combinations satisfying the constraint 
A*2 > A*i — f^s which, though being affected by larger statistical errors, are safe from the 
effects of the potentially divergent chiral logarithms. In table [3] we quote as a systematic 
error within each fit the uncertainty associated with the determination of the lattice spacing 
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Fit 


Polynomial 


PQChPT 


C-PQChPT 


<F (MeV) 
mfS (MeV) 
m s /m ud 
fx (MeV) 


4.07(9)(33) 
109(2)(9) 
26.7(2)(0) 
158.7(11)(89) 
1.214(8)(0) 


3.82(15)(25) 
107(3)(7) 
27.9(2)(0) 
160.2(15)(54) 
1.225(H)(0) 


3.74(13)(21) 
102(3)(6) 
27.4(3) (0) 
161.8(10)(0) 
1.238(7)(0) 



Table 3: Results for the light quark masses and pseudoscalar decay constants, in physical 
units, as obtained from the polynomial, PQChPT and C-PQChPT fits respectively, by 
analysing only the combinations of quark masses satisfying the constraint [i 2 > \i\ = Us- 
The quoted errors are statistical (first) and systematic (second), the latter coming from the 
uncertainties in the determination of the lattice scale and of the quark mass renormalization 
constant. 

and of the quark mass renormalization constant. 

In order to derive our final estimates for the quark masses and decay constants, we 
perform a weighted average of the results of the three analyses presented in table E] and 
conservatively add the whole spread among these results to the systematic uncertainty. In 
this way, we obtain as our final estimates of the light quark masses the results 

mj? (2 GeV) = 3.85 ± 0.12 ± 0.40 MeV , mf E (2 GeV) = 105 ± 3 ± 9 MeV , (16) 

and the ratio 

ms/mud = 27.3 ± 0.3 ± 1.2 , (17) 

where the first error is statistical and the second systematic. 

For the kaon decay constant and the ratio fx / f n we obtain the accurate determinations 

f K — 161.7 ± 1.2 ± 3.1 MeV , f K /U = 1-227 ± 0.009 ± 0.024 . (18) 

It is interesting to compare our result for the strange quark mass with other lattice 
QCD determinations of the same quantity. This comparison is illustrated in fig. HI 

The HPQCD-MILC-UKQCD Collaboration, using the MILC extensive simulations of 
lattice QCD performed with Nf = 2 + 1 dynamical improved staggered fermions, initially 
quoted the result m^ s (2 GeV) = 76(3) (7) MeV [H], significantly lower than our prediction 
in eq. (fl6|) . In |14j . the quark mass renormalization constant was determined using one- 
loop perturbation theory. The two-loop calculation has then led to a significant increase of 
the quark mass estimate [21] , and the most recent determination presented by MILC now 
reads m ] f s (2 GeV) = 90 (5) (4) MeV [25]. Recently, a similar result has been also obtained 
by the CP-PACS and JLQCD Collaborations, using (9(a)-improved Wilson fermions with 
Nf = 2 + 1 and implementing the quark mass renormalization at one loop: m^ ls (2 GeV) = 
91.1 (jlg 4 2 6 ) [26]. It is worth noting that this result is perfectly consistent with the previous 
Nf = 2 determinations of the same quantity obtained by the two collaborations [T7J ITS] . 
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Figure 4: Lattice QCD determinations of the strange quark mass obtained from simulations with 
2 JTTJ-ISf and N f = 2 + 1 JH El H3 \M dynamical fermions. The PDG average (from 
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f 



lattice only) [27} is also shown for comparison. 



In the present analysis, we find that the use of non-perturbative renormalization plays 
a crucial role in the determination of the quark masses. The estimate Zp l ~ MOM (l/a) = 
0.39(1) (2) obtained with the RI-MOM method is in fact significantly smaller than the 
prediction Zp FT (l/a) ~ 0.57(5) given by one-loop boosted perturbation theory (in the 
same RI-MOM renormalization scheme) |15j . Had we used the perturbative estimate of 
Z P we would have obtained mjjf (2 GeV) = 2.63 ± 0.08 ± 0.36 MeV and mf(2 GeV) = 
72 ± 2 ± 9 MeV. As shown in fig. HJ our prediction for the strange quark mass in eq. (TT6|) 
is in good agreement with other determinations based on a non-perturbative evaluation 
of the mass renormalization constant. These include the results obtained by ALPHA, 
mf (2 GeV) = 97(22) MeV [T3], by SPQ CD R, mf(2 GeV) = 101 (8) (tf) MeV [20], by 
QCDSF-UKQCD, mf(2 GeV) = 119(5)(8) MeV from the vector Ward identity [21] and 
mf s (2 GeV) = 111(6)(4)(6) MeV from the axial one [22], and by RBC, mf s (2 GeV) = 
119.5(56) (74) MeV [23]. It is often found that, in lattice determinations of quark masses, 
implementing a non-perturbative renormalization method has an impact that can be larger 
even than the quenching effect. We believe that this observation should be always kept in 
mind, particularly when the lattice results for quark masses are considered for producing 
final averages. 

Our result for the ratio fx/fn is compared in fig. [5] with other recent lattice determi- 
nations based on simulations with Nf = 2 and Nf = 2 + 1 dynamical fermions. 

Our calculation and those based on the MILC improved staggered gauge configura- 
tions are the only ones in which light quark masses significantly lower than m s /3 have 
been simulated (m q ~ m s /6 for our lightest quark mass). Therefore, it is interesting 
to compare our determination of fx/fn with the more recent results quoted by MILC, 
f K /f n = 1.208(2)(t[ 4 ) [25], and by HPQCD-UKQCD, f K /f n = 1.189(7) [33J. Despite the 
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Figure 5: Lattice QCD determinations of the ratio fx/fn obtained from simulations with Nt = 
2 [ig HHF and N f = 2 + 1 {25f, j2^-j3^ dynamical fermions. The results are also compared 
with the PDG 2006 average \27] and with the average based on the updated determination ofV us 
from Ku, decays \34\j . 

strange quark is still quenched in our simulation, and our results are still obtained at a 
single value of the lattice spacing, we find the agreement between these determinations 
quite satisfactory. In order to better quantify the size of discretization effects, which are 
of 0(a 2 ) in the present calculation, we plan to extend the simulation to other two values 
of the lattice spacing (corresponding to (3 — 3.8 and (3 = 4.05). This should also allow us 
to eventually perform the extrapolation to the continuum limit. 

Our result for the ratio fx/ fir can be combined with the experimental measurement of 
F(K — > nu^)) /Y{j[ — > jiu^ij)) [27] to get a determination of the ratio |V^ s |/|V^d| [35] . 
We obtain 

\V us \/\V ud \ = 0.2251(5)(47), (19) 

where the first error is the experimental one and the second is the theory error com- 
ing from the uncertainty on f K j f w . Eq fTl9|) . combined with the determination \V u d\ = 
0.97377(27) [36J from nuclear beta decays, yields the estimate 

\V US \ = 0.2192(5)(45), (20) 

in agreement with the value extracted from Kg 3 decays, |V^ S | = 0.2255(19) [34] , and leads 
to the constraint due to the unitarity of the CKM matrix 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 - 1 = (-3.7 ± 2.0) ■ 10~ 3 . (21) 

Acknowledgements 

We thank D. Becirevic, G. Martinelli and G.C. Rossi for useful comments and discussions. 



15 



The computer time for this project was made available to us by the John von Neumann- 
Institute for Computing on the JUMP and Jubl systems in Jiilich and apeNEXT system in 
Zeuthen, by UKQCD on the QCDOC machine at Edinburgh, by INFN on the apeNEXT 
systems in Rome, by BSC on MareNostrum in Barcelona (www.bsc.es) and by the Leibniz 
Computer centre in Munich on the Altix system. We thank these computer centres and 
their staff for all technical advice and help. 

This work has been supported in part by the DFG Sonderforschungsbereich/Transregio 
SFB/TR9-03, the EU Integrated Infrastructure Initiative Hadron Physics (I3HP) un- 
der contract RII3-CT-2004-506078 and the EU contract MRTN-CT-2006-035482, "FLA- 
VIAnet". We also thank the DEISA Consortium (co-funded by the EU, FP6 project 
508830), for support within the DEISA Extreme Computing Initiative (www.deisa.org). 
R.F., V.L. and S.S. thank MIUR (Italy) for partial financial support under the contracts 
PRIN06. V.G. and D.P. thank MEC (Spain) for partial financial support under grant 
FPA2005-00711. M.P. acknowledges financial support by an EIF Marie Curie fellowship of 
the European Community's Sixth Framework Programme under contract number MEIF- 
CT-2006-040458. 

Appendix 

In this appendix we collect in tables H] and [5] the values of the pseudoscalar meson masses 
and decay constants obtained at the various combinations of simulated sea and valence 
quark masses. 
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Table 4: Values of the pseudoscalar meson masses aM PS (fi s , /ii, /^) for the various combi- 
nations of simulated sea and valence quark masses. 
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0.0150 


0.0771 


V 4 ) 


0.0771 


[A) 


0.0782(5) 


0.0780 


[A) 


0.0789 


'5) 


0.0085 


0.0220 


0.0802 


V 4 ) 


0.0804 


'A) 


0.0814(5) 


0.0812 


[A) 


0.0821 


'5) 


0.0085 


0.0270 


0.0822 


;4) 


0.0825 


[A) 


0.0834(5) 


0.0832 


[A) 


0.0841 


'5) 


0.0085 


0.0320 


0.0841 


V 4 ) 


0.0844 


[A) 


0.0853(5) 


0.0850 


[A) 


0.0859 


;s) 


0.0100 


0.0100 


0.0753 


V 4 ) 


0.0753 


[A) 


0.0764(5) 


0.0763 


[A) 


0.0771 


;s) 


0.0100 


0.0150 


0.0779 


V 4 ) 


0.0780 


[A) 


0.0791(4) 


0.0789 


[A) 


0.0799 


;s) 


0.0100 


0.0220 


0.0811 


V 4 ) 


0.0813 


[A) 


0.0823(4) 


0.0820 


[A) 


0.0830 


;s) 


0.0100 


0.0270 


0.0831 


V 4 ) 


0.0834 


[A) 


0.0843(4) 


0.0840 


[A) 


0.0850 


'5) 


0.0100 


0.0320 


0.0850 


V 4 ) 


0.0853 


[A) 


0.0862(5) 


0.0859 


[A) 


0.0868 


'5) 


0.0150 


0.0150 


0.0806 


V 4 ) 


0.0808 


[A) 


0.0817(4) 


0.0814 


[A) 


0.0825 


V 5 ) 


0.0150 


0.0220 


0.0838 


V 4 ) 


0.0841 


[A) 


0.0849(4) 


0.0846 


[A) 


0.0857 


;s) 


0.0150 


0.0270 


0.0858 


;4) 


0.0862 


[A) 


0.0870(4) 


0.0866 


[A) 


0.0877 


'5) 


0.0150 


0.0320 


0.0877 


;4) 


0.0881 


[A) 


0.0889(4) 


0.0885 


[A) 


0.0896 


;s) 
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Table 5: Values of the pseudoscalar decay constants afps(i^s, A*i, ^2) f° r ^ e various com- 
binations of simulated sea and valence quark masses. 



[2] 
[3] 

[4] 

[5] 



S. R. Sharpe and J. M. S. Wu, Phys. Rev. D 71 (2005) 074501 arXiv:hep-lat/0411021 



References 

[1] Ph. Boucaud et al [ETM Collaboration], Phys. Lett. B 650 (2007) 304 
[arXiv:hep-lat/07010l2] . 

R. Frezzotti and G. C. Rossi, JHEP 0410 (2004) 070 [ arXiv:hep-lat /0407002| . 

A. M. Abdel- Rehim, R. Lewis, R. M. Woloshyn and J. M. S. Wu, Phys. Rev. D 74 
(2006) 014507 |arXiv:hep-lat/0601036| . 

M. Foster and C. Mich ael [UKQCD Collaboration], Phys. Rev. D 59 (1999) 074503 
|arXiv:hep-lat/981002l] . 

C. McNeile and C. Mich ael [UKQCD Collaboration], Phys. Rev. D 73 (2006) 074506 
|arXiv:hep-lat /0603007] . 

ETM Collaboration, in preparation. See also C. Urbach, plenary talk given at LAT- 
TICE 2007, |http : //www . physik . uni-regensburg . de/lat07~j 

[7] R. Frezzotti and G. C. Rossi, JHEP 0408 (2004) 007 [ arXiv:hep-lat/0306014| . 

[9 
[10 

[11 

[12 
[13 
[14 

[15 

[16 
[17 



R. Frezzotti, G. Martin elli, M. Papinutto and G. C. Rossi, JHEP 0604 (2006) 038 
|arXiv:hep-lat /0503034] . 

L. De l Debbio, L. Giusti, M . Luscher, R. Petronzio and N. Tantalo, JHEP 0702 (2007) 
082 [arXiv:hep-lat/0701009| . 

S. R. Sharpe, Phys. R ev. D 56, 7052 (1997) [Erratum-ibid. D 62, 099901 (2000)] 
|arXiv:hep-lat/9707018] . 



J. Bijnens and T. A. Lahde, Phys. Rev. D 72 (2005) 074502 [arXiv:hep-l at/0506004|. 

D. Becirevic and G. Villadoro, Phys. Rev. D 69 (2004) 054010 | arXiv:hep-lat/0311028] . 

C. Aubin et al. [HP QCD Collaboration], Phys. Rev. D 70 (2004) 031504 
|arXiv:hep-lat/0405022] . 

ETM Collaboration, in preparation. See also P. Dimopoulos, talk given at LATTICE 
2007, http : //www . physik . uni-regensburg . de/lat07 , 

K. G. Chetyrkin and A. Retey, Nucl. Phys. B 583 (2000) 3 | arXiv:hep-ph/9910332] . 

A. Ali Khan et al. [CP-PA CS Collaboration], Phys. Rev. D 65 (2002) 054505 [Erratum- 
ibid. D 67 (2003) 059901] larXiv:hep-lat/0105015] . 



[18] S. Aoki et al. [JLQ CD Collaboration], Phys. Rev. D 68 (2003) 054502 
[arXiv:hep-lat/02 12039] . 



19 



[19] M. Delia Morte, R. Hoffmann, F. Knechtli, J. Rolf, R. Sommer, I. Wet- 
zo rke and U. W olff [ALPHA Collaboration], Nucl. Phys. B 729 (2005) 117 
|arXiv:hep-lat /05070"35] . 



[20 
[21 

[22 
[23 

[24 

[25 

[26 

[27; 

[28 
[29 

[30 

[31 
[32 

[33 

[34 

[35 
[36 



D. Becirevic et al, Nucl. Phys. B 734 (2006) 138 | arXiv:hep-lat/0510014| . 

M. Gockeler, R. Horsley, A. C. Irving, D. Pleiter, P. E. L. Rakow, G. Schier- 
holz and H. Stuben [QCDSF Collaboration], Phys. Lett. B 639 (2006) 307 
|arXiv:hep-ph/0409312| . 

M. Gockeler et al, Phys. Rev. D 73 (2006) 054508 |arXiv:hep-lat/0601004| . 

T. Blum, T. Doi, M. Hayakawa, T. Izubuchi and N. Yamada, larXiv:0708.0484l [hep- 
lat]. 

Q. Mason, H. D. Trottier, R. Horgan, C. T. H Davies and G. P. Lep age [HPQCD 
Collaboration], Phys. Rev. D 73 (2006) 114501 [arXiv:hep-ph/0511160] . 

C. Bernard et al [MILC Collaboration], PoS LAT2006 (2006) 163 
|arXiv:hep-lat /0609053] . 

T. Ishikawa et al. [JLQCD Collaboration], larXjyj07I)44M/] [hep-lat]. 
W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 



Y. Aoki et al, Phys. Rev. D 72 (2005) 114505 [arXiv:hep-lat/04110 06|. 

C. Aubin et al. [M ILC Collaboration], Phys. Rev. D 70 (2004) 114501 
|arXiv:hep-lat /0407028] . 

S. R. Be ane, P. F. Bedaque, K . Orginos and M. J. Savage, Phys. Rev. D 75 (2007) 
094501 |arXiv:hep-lat/0606023] . 

T. Ishikawa et al, PoS LAT2006 (2006) 181 | arXiv: hep-lat /06 10050] . 

C. Allton et al. [RBC a nd UKQCD Collaborations], Phys. Rev. D 76 (2007) 014504 
|arXiv:hep-lat/0701013] . 

E. Follana, C. T. H. Davies, G. P. Lepage and J. Shigemitsu [HPQCD Collaboration], 
larXiv:0706.l726l [hep-lat]. 



G. Isidori, conference summary talk at KAON'07, http://www.lnf.infn.it/ 



conf erence/kaon07; see also G. Isidori. larXiv:070 9.2438 [hep-ph]. 



W. J. Marciano, Phys. Rev. Lett. 93 (2004) 231803 |arXiv:hep-ph~70 402299 1 . 

W. J. Marciano an d A. Sirlin, Phys. Rev. Lett. 96 (2006) 032002 
|arXiv:hep-ph/05 10099] . 



20 



